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Описывается процедура идентификации параметров объекта управления, основанная на применении функции чувствительно- 
сти. Предлагается алгоритм локальной оптимизации, позволяющий идентифицировать неизвестные параметры из условия ми- 
нимизации квадратичного критерия невязки наблюдаемых переменных и оценки состояния на дискретных интервалах времени. 
Рассматривается пример алгоритма идентификации для определения динамического и статического моментов электропривода 
постоянного тока. 


Введение 

Жесткие требования к качественным показате- 
лям процесса движения многих технологических 
объектов определяют необходимость применения 
адаптивного управления исполнительными элек- 
тромеханическими системами. Для синтеза адап- 
тивной системы управления часто используется 
эталонная модель электромеханического объекта, 
которая предполагает достаточно точное знание 
его параметров [1]. При этом практический инте- 
рес представляет собой задача идентификации па- 
раметров электромеханического объекта, изме- 
няющихся детерминированным образом в преде- 
лах области ограниченных значений. В теории 
идентификации параметров детерминированных 
систем значительное место занимают методы иден- 
тификации, основанные на применении функции 
чувствительности [2] . Это объясняется следующи- 
ми свойствами методов: во-первых, они имеют до- 
статочно быструю сходимость процедур идентифи- 
кации; во-вторых, высокую точность результатов и 
самым определяющим фактором является их уни- 
версальность. Данные методы применимы как к 
линейным системам, так и к сложным нелинейным 
системам уравнений, описывающим поведение ре- 
ального объекта в пространстве состояний. 

Постановка задачи 

В следящих электромеханических системах воз- 
никает задача оценки, с последующей компенса- 
цией, влияния на характеристики регулируемого 


процесса таких параметров исполнительных при- 
водов, как динамический и статический моменты 
механической нагрузки. Одним из вариантов ре- 
шения такой задачи является включение в контур 
управления идентификатора параметров [1, 3, 4]. 
Кроме того, современные цифровые системы упра- 
вления электромеханическими объектами позво- 
ляют расширить математическое обеспечение и ис- 
пользовать более эффективные алгоритмы иденти- 
фикации неконтролируемых параметров в режиме 
реального времени. 

Рассмотрим алгоритм идентификации на примере 
электромеханического объекта - электропривода по- 
стоянного тока с упругой механической связью, пара- 
метры которой изменяются неизвестным образом в 
определенных пределах. Динамика электропривода с 
полупроводниковым преобразователем и двухмассо- 
вым механическим звеном может быть описана систе- 
мой дифференциальных уравнений [3, 4]: 

(Т п р + \)е п =к г к п іі у , 


а ^=м-м 12 -м л , 

Асо 2 _ 

°2 ^ 1Ѵ1 12 ІУ1 с29 

М\2 — С] 2 {Фі — Ф 2 )і 


( 1 ) 


где ку, к„ - коэффициенты передачи усилителя и 
преобразователя; Т„ - электромагнитная постоян- 
ная времени преобразователя; е„ - выходное на- 
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пряжение преобразователя; и у - управляющее воз- 
действие (выходной сигнал регулятора); Т я - элек- 
тромагнитная постоянная времени якорной цепи; 
М— электромагнитный момент двигателя; К„ — ак- 
тивное сопротивление якоря двигателя; к д - кон- 
структивный коэффициент двигателя; І ь / 2 - мо- 
менты инерции двигателя и механической нагруз- 
ки; со ь со 2 - угловые частоты вращения вала двига- 
теля и механической нагрузки; М а - упругий мо- 
мент кинематической цепи; М А , М й - момент со- 
противления на валу двигателя и механической на- 
грузки; с 12 - коэффициент упругости кинематиче- 
ской цепи; ф ь ф г — углы поворота вала двигателя и 
механической нагрузки. 

Для данного объекта в качестве переменных со- 
стояния принимаются: напряжение е п , моменты М , 
М А , М а . Предполагается, что измеряемыми пере- 
менными является полный вектор состояния, а 
входными воздействиями - выходной сигнал и у ре- 
гулятора, моменты М А и М а . Вследствие нестабиль- 
ности механической нагрузки необходимо иденти- 
фицировать значения следующих параметров и пе- 
ременных: моменты инерции ./, и ^ ъ моменты со- 
противления М А и М а . 

Систему ур. (1) приведём к нормальному виду с 
учетом следующих обозначений: 
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В результате запишем: 
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Таким образом, рассматриваемый электропри- 
вод, имеющий несколько входов и выходов, явля- 
ется многомерной линейной системой (2). Для 
упрощения его математического описания удобно 
воспользоваться моделью в векторно-матричной 
форме. При этом п переменных состояния заменя- 


ются одной обобщённой переменной — вектором 
состояния X, а г входных переменных — вектором 
входов II. Тогда система уравнений (2) представля- 
ется уравнением в векторно-матричной форме 


Х(і) = А(Ѳ)Х(() + В(Ѳ)Щ(), і е [і„Т], (3) 
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Ѳ - вектор неизвестных параметров: 


Ѳ 1 = Р ѳ 2 =а~, ѳ 3 =м с1 , ѳ 4 = м с1 . 


Синтез алгоритма идентификации 

Задача идентификации неизвестных параметров 
заключается в следующем: по известным наблюде- 
ниям 2(і)=Х(і), ?б[/ 0 , 7] определить оценку вектора 
значений Ѳ, минимизируя критерий /(©)—> тіп ка- 
чества в пространстве параметров Ѳ. Решение по- 
добных задач рассматривалось в работах [5-7]. 

В случае полной информации с выхода объекта 
2{і)=Х(і), ?е[/ 0 , 7] критерий качества идентифика- 
ции параметров выбирается в виде: 

У(Ѳ)= \(2(і)-Х(і,®)) Т {2(і)-Х(і,&))сІі, 

і (4) 

& = (Ѳ Ѵ ...,Ѳ 4 ) Т . 

При подстановке решения ур. (3) в (4), крите- 
рий качества является нелинейной функцией от 
вектора неизвестных параметров Ѳ, и минимиза- 
ция критерия (4) приведёт к системе нелинейных 
уравнений относительно Ѳ. Очевидно, что реше- 
ние системы нелинейных уравнений будет пред- 
ставлять значительные трудности. Для решения за- 
дачи целесообразно представить Х(1,Ѳ) в окрестно- 
сти равновесного состояния Х(і.Ѳ"), которое явля- 
ется решением уравнения 

х(і 0 ,ё°) = х(ір. 

Ограничимся двумя членами разложения 

Я у „ 

ЛГ(/,Ѳ) « ЛГ(/,Ѳ°)+ (Ѳ-Ѳ 0 ), (5) 

5© ѳ= в» 

где Ѳ'-Ѳ -АѲ 1 '; Ѳ° - нулевое приближение; Ѳ 1 - 
оценка вектора неизвестных параметров, получен- 
ная на первой итерации; частная производная 
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называется функцией чувствительно- 
го® ѳ=ѳ* 

сти параметров и обозначается как а>(1,&), где 
к - номер итерации процесса идентификации. При 
этом переменные Д/,Ѳ°), ет(?,Ѳ°) удовлетворяют 
следующим уравнениям (при к= 0): 

Х(1,®) = А(®)Х(і,®) + В(®)Ѵ(і), 

Х(І к ,®) = Х к , 

= Х((,®) + 

д ® ѳ=ѳ‘ 

+ДѲ‘Ж?,Ѳ‘)+ а ^ Щі), (6) 

0 ^ ѳ=ѳ* 

®(4,©‘) = о. 

При подстановке (5) в (4) критерий качества 
принимает вид 

ДѲ°) = | (2(1) - Х(і, Ѳ' 1 ) - со(і, Ѳ°)ЛѲ' , ) Т х 

1 0 

*(2(1)-Х(1,®)-со(1,®) АѲ°) Л. (7) 

Необходимое условие минимума функции (7) 
по АѲ° примет вид 

АѲ° = 

1 0 

= | со((,Ѳ°) т (2 — Х((,Ѳ°)) ск. 

1 0 

В результате решения данной алгебраической 
системы определяется приращение АѲ°. 

На основе вышеприведённых формул можно 
построить следующую итерационную процедуру. 
Пусть на А;-ой итерации получено значение оценки 
вектора параметров &. Тогда решение Х(1,Ѳ) мож- 
но разложить в окрестности Х((,&): 

Х(і, ®)*Х(і,Ѳ к ) + ю ((, Ѳ к )((, Ѳ к )АѲ а , (8) 

где Д/,Ѳ‘), со(1,&) удовлетворяют ур. (6). 

Подставив (8) в критерий (4) и, минимизируя его 
по А&, получим следующую алгебраическую систему 

\(о(і,ё к ) т (о(і,ё к )сіі- аѳ * = 

‘о 

-\со т (1,® к )(2-Х(1,®)) ск. 

1 0 

Обозначая 

Ф = |й) г (?,Ѳ')юН,Ѳ‘)Л,Т / = 

'о 

= ]со т (і,®)(2-Х(і,®))ск, (9) 

1 0 


можно записать 

ФДѲ А =Ф. (10) 

В результате решения системы (10), можно вы- 
числить новое (ДН)-ое приближение идентифи- 
цируемых параметров по формуле 

Ѳ а+1 =Ѳ а +АѲ а . (11) 

Пусть выполнено к итераций и получена оцен- 
ка вектора параметров на к- ом шаге. Тогда на осно- 
ве вышеприведённых формул можно записать сле- 
дующий алгоритм: 

1. Используя наблюдения 2(1), формируются ма- 
трицы Ф и Ф по формулам (9). 

2. Решается алгебраическая система (10) для на- 
хождения приращения АѲ*. 

3 . Вычисляется следующая (к+ 1 ) -ая оценка вектора 
идентифицируемых параметров по формуле (11). 
Рассмотрим модификацию описанного метода 

идентификации на основе функции чувствитель- 
ности, называемую методом локальной оптимиза- 
ции. Для того чтобы не запоминать всю реализа- 
цию 2(1), можно минимизировать критерий каче- 
ства (4) последовательно во времени. С этой целью 
на интервале времени [4,7] построим равномер- 

Т — і 0 

ную сетку с узлами і к =і 0 +к — — — ,к = О, А. 
Критерий качества (4) представим в виде суммы 

к = ХЛ, г Д е Л = } \2(1)-Х(1,®)^ (к. 

При этом последовательность векторов оценок 
Ѳ°, Ѳ 1 , Ѳ ѵ формируется следующим образом. Век- 
тор оценок ѲД1 находится по формуле 

Ѳ Ж =Ѳ 4 +ЛѲ Ж , к= 0,7Ѵ — 1 , Ѳ° - нулевое приближе- 
ние. Приращение АѲ‘ определяем из условия ми- 
нимума критерия А к -> тіп . Подставим в критерий 

АѲ к 

./, значения вектора Х(!,Ѳ) в виде 

Х(і,Ѳ)*Х((,Ѳ к )+ дХ ( 1,Ѳ) ■ ДѲ А , 

гѲ ѳ=( §* 

Л7-/ Р\к\ / Р\к \ 

где Х(і,Ѳ )и (о(і, Ѳ ) = — — удовлетворяют 
ур. (6). 

Минимизируя критерий по значениям ДѲ*, 
получим: 

| ю(і,Ѳ к Усо(і,®)ск- АѲ‘ = 

І к 

= I о/ (1,®)(2-Х (!,&)) Ас. 

І к 

Из данной алгебраической системы определя- 
ется приращение ДѲ*. 

Введём следующие обозначения: 
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а) 



Ь) 


Рисунок. Изменение значений оценок параметров Ѳі Ѳ 2 : а) с 
увеличением числа итераций и Ь) с течением вре- 
мени 
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2=| со т (1,&)со(1,&)сІІ, 

І к 

0 = | ( 12 ) 

І к 

Пусть выполнено к итераций и получена оцен- 
ка вектора АѲ* параметров на к - ом шаге. С учётом 
обозначения 

2АѲ* = О., (13) 

алгоритм идентификации будет заключаться в сле- 
дующем. 

1. Используя наблюдения Д і), формируем матри- 
цы 2 и О по формулам (12). 

2. Решается алгебраическая система (13) для на- 
хождения приращения ДѲ\ 

3. Вычисляется следующая (к+ 1)-ая оценка векто- 
ра параметров по формуле (11). 

Результаты моделирования 

В качестве примера моделирования рассмотрим 
алгоритм локальной оптимизации. На рисунке при- 
ведены результаты идентификации параметров Ѳ ь 
Ѳ>. Известны истинные значения параметров: 
Ѳ ““=454, 545, Ѳ“"=263,157 (т.к. /Д=0,0022 и 
/“"=0,0038 соответственно). В результате работы ал- 
горитма получены следующие значения: Ѳ,=454,449, 
02=263,8 (/[=0,0022 и /=0,00379 соответственно). 

Выводы 

Анализ результатов моделирования отражает 
асимптотический характер изменения значений оце- 
ниваемых параметров. Сходимость параметров к ис- 
тинным значениям осуществляется за сравнительно 
короткое время. Разработанный алгоритм можно ис- 
пользовать при построении адаптивных систем упра- 
вления сложными электромеханическими объектами. 
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